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Abstract 

A connection between Romanovski polynomials and those polyno- 
mials that solve the one-dimensional Schrodinger equation with the 
trigonometric Rosen-Morse and hyperbolic Scarf potential is estab- 
lished. The map is constructed by reworking the Rodrigues formula 
in an elementary and natural way. The generating function is summed 
in closed form from which recursion relations and addition theorems 
follow. Relations to some classical polynomials are also given. 
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1 Introduction and review of basic results 

Romanovski polynomials were discovered in 1884 by Routh [1] in the form of 
complexified Jacobi polynomials on the unit circle in the complex plane and 
were then rediscovered as real polynomials by Romanovski [2] in a statistics 
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framework. Recently real polynomial solutions of the Scarf [3] and Rosen- 
Morse potentials [I] in (supersymmetric) quantum mechanics were recog- 
nized [5] to be related to the Romanovski polynomials. 

Here we apply to Romanovski polynomials a recently introduced natural 
method of reworking the Rodrigues formula |6| that leads to connections with 
other polynomials. 

The paper is organized as follows. In the Introduction we define the com- 
plementary polynomials Q[^~ a \x) then establish both a recursive differential 
equation satisfied by them and the procedure for the systematic construction 
of the Q^~ a \x), derive their Sturm-Liouville differential equation (ODE), 
their generating function and its consequences, all based on the results of 
ref. [5J. Section 2 deals with a parameter addition theorem, Sect. 3 with 
orthogonality integrals and Sect. 4 with connections of these Romanovski 
polynomials to the classical polynomials. Sect. 5 deals with further applica- 
tions using auxiliary polynomials. 

Definition. The Rodrigues formula that generates the polynomials is 
given by 

pM/ \ 1 d * r m nil °j x )l dlw o( x ) , n i fi\ 
p i {x) = —^ — [wi(x)a(xy] = — 7T , 1 = 0,1,.... 1 

Wi(x) dx 1 Wq(x) dx 1 

where cr(x) = 1 +x 2 is the coefficient of y" of the hypergeometric ODE (1) of 
ref. [6J that the polynomials satisfy. The variable x is real and ranges from 
— oo to +oo. The corresponding weight functions 

Wl ( x ) = a{x)-^ l+a+ ^e- acot ~ lx = a(x)~ l w (x) (2) 

depend on the parameters a, a that are independent of the degree I of the 
polynomial Pi' a \x). o 

Lemma. The weight functions of the P} (x) polynomials satisfy Pear- 
son's first-order ODE 

a(x) =[a- 2x(l + a + l)]w t (x). (3) 

Proof. This is straightforward to check using dcot~ 1 x/dx = —l/a(x). o 

We now apply the method of ref. [6J and introduce the complementary 
polynomials Q ( ^ l ~ a \x) defining them inductively according to the Rodrigues 
representation (see Eq. (5) of ref. [6]) 

p^ a \x) = J f^M x y~ v Mx)Q { °'~ a) (x)l ^ = 0,1,...,/. o (4) 
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For Eq. PJ to agree with the Rodigues formula (PQ) for v = requires 
Qo (x) = 1. Then comparing Eq. ([3]) with Eq. (4) of ref. [5j we find 
the coefficient t(x) = a — o~'(x)(a + I) of y' in the ODE of the polynomi- 
als. Comparing instead with the ODE (37) of ref. [5] gives their parameter 
(3 = —a. 

We now identify the polynomials of ref. JBj with those defined in Eq. O) 

V u (x;l)=Q^- a \x), \ >v. (5) 

We will show below that the polynomials Q^~ a \x) are independent of the 
parameter /. o 

Theorem 1.1. Q^~ a \x) is a polynomial of degree v that satisfies the 
recursive differential relation 

Q<%- a \x) = a(x) ^ g ^J )(x) + [t(x) + 2x(l -v- \)\Q^~ a \x) 

= a{x) W » ^ 1 ! +[a-2x(a + u + l)]Q^- a \x), 

v = 0,1,.... (6) 

Proof. The inductive proof of Theorem 2.2 of ref. [6] applied to the poly- 
nomial Q^~ a \x) proves this theorem, and Eq. agrees with Eq. (76) of 
ref. [5] provided their parameter (3 = — a in ref. [6]. Since Eq. ([6]) is indepen- 
dent of the parameter I, so are the polynomials Q ( ^ l ~ a \x) that are generated 
from it. o 

Comparing the recursive ODE ([6]) with one stated in [2] leads us to the 
identification of our polynomials 

Qi a ' k - m \x) = Mm,x), Q^- a \x) = Ma + v,x) (7) 

as a Romanovski polynomial (with its parameter depending on its degree), 
and comparing with Eq. (69) of ref. [5], 

Qi a >- a \x) = Ri a ^\x), (3 = -a.o (8) 

Notice that the parameter a {—u in ref. [2]) is suppressed in Romanovski's 
notation. The fact that the integer index v of the complementary polynomi- 
als occurs in the parameter m of the Romanovski polynomials is sometimes 
disadvantageous (for orthogonality), but also occasionally a definite advan- 
tage (for the generating function). Moreover, 

V l (x-l) = Q^- a \x) = K l Cl a '- a \x), (9) 
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where K\ is a normalization constant and the Cf"' - (x), after a change of 
variables, become part of the solutions of the Scarf and Rosen-Morse poten- 
tials in the Schrodinger equation [3], [I]. However, for C\ a ~ a \x) to be part of 
the solution of the Schrodinger equation with the trigonometric Rosen- Morse 
potential requires a = a\ = ^ [3]; but the polynomials may be defined for 
the general parameter a. 

Recursion relations and recursive ODEs are practical tools to systemati- 
cally generate the polynomials. 

Theorem 1.2. The polynomial Q^~ a \x) satisfies the basic recursive ODE 

dQ \ {X) = -*(2a + v + l)Q { ^ a \x) = -KQtr\x). (10) 

Proof. Eq. f TTUT) follows from a comparison of the recursive ODE ([6]) with 
a three-term recursion relation as outlined in Corollary 4.2 of ref. [6], is 
ODE (32) of ref. [6], and agrees with Eq. (75) of ref. [5] provided their 
(3 = —a, which is consistent with our previous statements, o 

Thus, taking a derivative of Qs*'~ a '(x) just lowers its degree (and index) 
by unity, up to a constant factor, a property the Romanovski polynomials 
share with all classical polynomials. 

Theorem 1.3. The polynomials Qi?~ a '{x) satisfy the differential equation 
of Sturm-Liouville type 

a(x f ^ 1 ] + [a - a'(x)(a + ") f* v ^ K) = -\ v Q^~ a) {x). (11) 
Proof. Substituting the basic ODE (flOl) in the recursive ODE ([6]) yields 

^(a-a), x a ( X ) ^Qu+l ( x ) 1 r / . . ,\ n &Qv+\ \ x ) / 10 \ 

Ql+i '{x) = — — [a- (a + u + 1)<t] (12) 

a v +\ ax ^u+i ux 

which, for v —>■ v — 1 is the ODE of the theorem. Again, the ODE is 
independent of the parameter I in Eq. (jSJ). o 

Theorem The polynomial Pi(x) satisfies the ODE 

d 2 P t (x) ,dPi(x) 

, 2 +r{x)— = -X l P l {x). (13) 

Proof. For v = I, we use Pi(x) = Vi(x; I) in the notation of ref. [6] to rewrite 
the recursive ODE (J6]) of ref. [6J and Eq. (flQj) as 

Vi(x;l) = a{x)V' l - l {x;l)+T{x)V l - 1 {x;l)=P l {x) 

= - a -fpl\x)- T -f P ^ (14) 
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which is the ODE (1) in ref. [5J for the polynomial Pi(x). (Note that t(x) is 
given right after Eq. (J3]).) o 

Theorem 1.5. The polynomials Q^~ a \x) satisfy the generalized Ro- 
drigues formulas 

g(«.- )( X ) = wr 1 (x)a(x)^^ 7 (a(x)^(x)Q!r'- a) (^)), 

// = 0,1,...,!/. (16) 

Proof. These Rodrigues formulas are those of Theorem 2.3 of ref. [Bj; they 
agree with Eqs. (72) and (73) of ref. [5] provided their (3 = —a, as we found 
earlier. Note that, from Eq. (j2j) the product wi(x)a(x) 1 does not depend on 
I, so there is no I dependence in Eqs. (115)161) . o 

The Q\° i ~ a ' (x) polynomials generalize the P n (x) in the sense of allowing 
any power of a(x) in the Rodrigues formula, not just a(x) n as for the P n (x). In 
other words, the Q^~ a \x) are associated P n {x) (or C^~ a \x)) polynomials, 
as in the relationship between Laguerre (Legendre) and associated Laguerre 
(Legendre) polynomials. 

A generalization of Eq. (fl5l) . 



Q<°.— 0(x) = ^T-TT Uxy l w»{x)) , I = 0, ±1, . . . , (17) 
w (x) dx u v ' 

just reproduces the same polynomial with a shifted parameter a — > a + I. 
The generating function for the Q ( ^ l ~ a \x) polynomials is defined as 

Q(x,y;a,-a) = ^ (x). o (18) 

The generating function is our main tool for deriving recursion relations. 
Theorem 1.6 The generating function can be summed in the closed form 

Wl (x)Q(x,y;a-a) = a{xf- 1 £ ^ML^L. ( <7(a . ) -(-+D e «cat^ - 

u=Q V - aX 

Q(x,y;a,-a) = (1 + x 2 ) a+ V cot_1 x [l + (x + y{l + x 2 )) 2 ]- (a+1) 

_ e -acot- 1 (x+j/(l+a; 2 )) ^ 
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[wi(x)a(x) l ~ 1 Q(x, y; a, -a) 



dy 



i 1 



m,f £ ^^Wxl-^e-"'"^"'""'^)). (20) 



, : _ /( [v •//)! d.r- " 

d*Q(x,ujc*-a) = Wi{x) -i a{xY -i [{l + {x + ya{x) y r ^) 

■ e- acot ~ 1{x+ya{x))2 Q^- a \{x + ya{x)f). (21) 

Both Taylor series converge if x and x + yo~(x) are regular points of the weight 
function. 

Proof. The first relation is derived in Theorem 3.2 of ref. [6] by substi- 
tuting the Rodrigues formula (j!5p in the defining series (fT8l) of the generating 
function and recognizing it as a Taylor series. The other follows similarly, o 

Theorem 1.7. The generating function satisfies the partial differential 
equation (PDE) 

dQ(x,y;a,-a) _ a(x)Q(x,y;a,-a) {a _ a) 

dy " 1 + [x + ya{x)f Ql ^ X + y ^ X >>- ™ 

Proof. This PDE is derived by straightforward differentiation in Theo- 
rem 3.3 of ref. [6] in preparation for recursion relations by translating the 
case \i — 1 in Eq. ([211) into a partial differential equation (PDE). o 

One of the main consequences of Theorem 1.7 is a general recursion rela- 
tion. 

Theorem 1.8. The Q)° i ~ a \x) polynomials satisfy the three-term recursion 
relation 

Q { Xi a \x) = [a - 2x{a + u + l)]Q^~ a \x) - ua(x)(2a + u + l)Q { ^ a \x). 

(23) 

Proof. Equation (j22|) translates into this recursion relation by substitut- 
ing Eq. (|T8|) defining the generating function and thus rewriting this as 

oo „.v—X 



;i + y 2 a\x) + 2xy) £ -1—Q<*-4( X ) 

v =i \y L )- 

- 2x(a + 1) - 2y(a + l)a(x)} £ ^jQ^~ a \x). o (24) 



" t^o v 
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Just like the recursive ODE (J2]), this recursion allows for a systematic 
construction of the Romanovski polynomials, in contrast to the Rodrigues 
formulas which become impractical for large values of the degree v. 

Theorem 1.9. The polynomials Q*£*'~ a '(x) satisfy the differential equation 
of Sturm- Liouville type 

^(^) W ^i(^) dQ ^. a)(a?) ) = -Ko-{x) l -»{x) Wl {x)Q^- a \x)- 

A„ = u(2a + u + l) . (25) 

Proof. This ODE is equivalent to the ODE ( TlTT) and agrees with Eq. (78) 
of ref. [5] if f3 = —a there. Note that the inductive proof in Theorem 5.1 in 
ref. [6] is much lengthier than our proof of Eqs. ( fmfT3i) . o 



2 Parameter Addition 

The multiplicative structure of the generating function of Eq. f[T9l involving 
the two parameters in the exponents of two separate functions, as displayed 
in 

allows for the following theorems. 

Theorem 2.1. The Q^~ a \x) polynomials satisfy the parameter addition 
relation 

^ 1 =o Vz/ i / 

Proof. This formula follows from the Taylor expansion 

E y —Ql?'- ai) (x)Qi7'~ a2) (x) = E y -^Q { N 1+a2 - ai - a2 ~ 1] ^)- (28) 

of the generating function identity 

Q(x, y; ai, -ai)Q(x, y; a 2 , -a 2 ) = Q(x, y; a x + a 2 , -ai - a 2 - 1). o (29) 
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Given the complexity of the polynomials, the elegance and simplicity of 
this relation are remarkable. 

Example. The case N = is trivial, and N — 1 becomes the additive 
identity 

g^ 1+a2 '- (ai+a2+1)) (x) = ai + a 2 - 2x( ai + a 2 + 2) 

= [ai - 2x( ai + 1)] + [a 2 - 2x(a 2 + 1)] 

= Q[ ai '~ ai \x) + Q[ a2 '- a2 \x). (30) 

The first case involving additive and multiplicative aspects of the polynomials 
is N = 2 which we decompose and multiply out as follows: 

1+a2 '~ ai ~ a2_1) = [ai + a 2 - 2x(ai + a 2 + 2)] [ ai + a 2 - 2x(ai + a 2 + 3)] 

- 2a(x)(a 1 + a 2 + 2) = {[«!- 2z(ai + 1)] 

+ [a 2 -2x(a 2 + l)]}{[a;i-2:r(ai + 2)] 

+ [a 2 - 2a;(a 2 + 1)]} - 2<r(a:)[(ai + 1) + (a 2 + 1)] 

= [ai - 2x(ai + l)][ai - 2x(ai + 2)] - 2a(x)(ai + 1) 

+ [a 2 - 2x(a 2 + 1)][« 2 - 2x(a 2 + 2)] - 2a(x)(a 2 + 1) 

+ [ai - 2x(a l + 1)] [a 2 - 2x[a 2 + 1)] 

+ [a 2 - 2x[a 2 + 1)] [ai - 2x(a 1 + 1)] 

= Q^ 1 '"" 1 ^ + Q^ a ' 2 ~ a ' 2 ' ) _|_ Q( a i'~ a i)Q( a 2,-a2) 

+ gS a2 '- a2) gS ai, ~ ai) . o (31) 

The addition theorem is consistent with the homogeneous polynomial theo- 
rem in the variables x, a, yfa (without using a = x 2 + l) which the polynomials 
satisfy and can be generalized to an arbitrary number of parameters. 

Theorem 2.2. The Qu^ a ~ a \x) polynomials satisfy the more general poly- 
nomial identity 

^(ai+o 2 H \-a n ,-(a 1 +a 2 -\ \-a n +n-l)) , \ 

e ^n* , - ,1 (4 (32) 

0<v j <N,v l +-+v n =N+n All U r j=l 

Proof. This follows similarly from the Taylor expansion of the product 
identity of n generating functions 

n 

II Q(x, y~, 0(j, —dj) = Q(x, y;ai + a 2 -\ h a n , -(a x + a 2 H h a n + n - 1)). 

(33) 
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As an application of the parameter addition theorem we now separate the 
two parameters a and a into two sets of simpler polynomials Q^~ a ^ and 
Q(a,i)_ enc ^ we eX p anc [ the generating functions in the identity 

Q(x, y; a, -a) = Q(x, y; 0, -a)Q(x, y; a, 1) (34) 

in terms of their defining polynomials. This yields 

Theorem 2.3. The Q^~ a \x) polynomials satisfy the decomposition iden- 
tity 

N / N\ 

Qp~ a \x) = E r)Q^- a Hx)Qpy v (x). (35) 

Proof. This identity follows from a Taylor expansion of the generating 
function identity (13~4"|) in terms of sums of products of polynomials involving 
only one parameter each, o 

Definition. The generating function 



a[cot irr-cot 1 (x+ya(x))] _ V Q(a,l) ^ (%Q^ 

u=0 V - 



defines the second set of the polynomials, while the first one will be treated 
in detail below upon expanding the polynomials Q^~ a \x) as finite sums of 
Gegenbauer polynomials in Sect. IV or finite power series in Eq. (|53|) . o 

We also note that Q^~ a \x) = K u C(°~ a \x), so the latter also have a 
Gegenbauer polynomial expansion. 

Theorem 2.4- The Q^~ a \x) polynomials satisfy the (parity) symmetry 
relation 

Ql~^)( x ) = (-iy Q ^-\- x ). (37) 

Proof. This relation derives from the generating function identity 

Q(-x,-y;-a,-a) = Q(x,y;a,-a) (38) 

which holds because acot" 1 (a; + ycr(x)) in the generating function Eq. fflBl 
stays invariant under a —>■ —a, x — > —x, y — > —y, and o~(—x) = a(x). o 
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3 Orthogonality Integrals 



This section deals with an application of the generating function to integrals 
that are relevant for studying the orthogonality of the polynomials. 

Definition. We define orthogonality integrals for the Q^~ a \x) polyno- 
mials by [5] 

Q^- a )( x )Q^'- a \x)e- acot ~ lx 
O a ; Q) = / dx-^ V /7T . = 0, 



W /._„ (J (r r .\{n+v)/2+a+2 ' ' 



r x \(ji+v)/2+a-\ 



> -3/2, n + v even, (39) 



while for \i + v odd there needs to be an extra ycr(x) in the numerator for 
the orthogonality integrals to vanish, o 

Thus, the Q^~ a \x) polynomials form two infinite subsets, each with gen- 
eral orthogonality, but polynomials from different subsets are not mutually 
orthogonal. While displaying infinite orthogonality, this property falls short 
of the general orthogonality of all classical polynomials. The Q^~ a \x) poly- 
nomials form a partition of the set of all Romanovski polynomials, as shown 
in Eq. (jSJ), with upper index dependent on the running index u, though. 
The Romanovski polynomials R^'^ix) with upper indices independent of 
the degree u, the running index, form another partition that has the finite 
orthogonality, as discussed in more detail in ref. [5]. The orthogonality of the 
C^° ra '(i) polynomials from the Schrodinger equation with a = a n as dis- 
cussed below Eq. is yet another form of orthogonality similar to that of 
hydrogenic wave functions, which also differs from the mathematical orthog- 
onality of associated Laguerre polynomials, the subject of Exercise 13.2.11 
in ref. [TJ. The orthogonality integrals of Eq. fl39l) suggest analyzing the 
following integral of the generating functions 



m = L^p^( g(x '^ ;a '- a) ) e " 



■a cot x 



y vi+»2 -oo Qi a '~ a \x)Q^- a \x)e 



-a cot 1 x 



V - / dx- 

,^=0 VlW. J-oo a ( x) (v 1+ u 2 )/2 + a + 2 



y - — o {a ' a) (40) 



OO yV!+U 2 
V1,V2=Q 1 Z 

which is written directly in terms of orthogonality integrals 0^~f> defined 
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in Eq. (I39p . On the other hand, we can express the integral as 







00 dx€~ aC °^ lx +2«cot 1 x— 2a cot 1 (x+y^/a) 
-00 a (a;)a+2[ 1 + y 2 + ^]2(a+l) 

00 d x gCicot -1 x-2ocot _1 (x+j/ v ^r) 
a a+2 [! + y 2 + ^]2(a+l) 



Z" 00 (ia; g-acot 1 x+2acot 1 (x—y^/a) 
+ X ^a+2 [! +y 2 _ ^]2(a+l) ' ( 41 ) 

which is manifestly not even in the variable y. If the Q^'^ polynomials 
were orthogonal, then the double sum in / of Eq. fj40l) would collapse to 
a single sum over normalization integrals multiplied by even powers of y, 
i.e. I would be an even function of y. This result shows that the Q^~ a ^ 
polynomials are not orthogonal in the conventional sense. In fact, the extra 
yfa in the orthogonality integrals O^'^L-i * s n °t built into the generating 
function. In other words, the fact that I(y) ^ I{—y) indirectly confirms 
that the Q^~ a ^ polynomials have more complicated orthogonality properties 
than the Romanovski polynomials with parameters that are independent of 
the degree of the polynomial, as discussed in more detail in ref . [5] . 

Let us next consider the special parameter a = and analyze similarly 
the integral 

dx y n s ™y u [°° dxQ^- a \x) 



Io{y) = L^xW> Q{x ^' A - a) = ^.L 



^ y v n (*>°) f°° dx 1 . f°° dx 



v\ u '° Jo a a+2 [1 + y 2 + Jo a a+2 [1 + y 2 - 2 m^ ' 



v =o w - l x ^ y ^ " u " l " ^ 

which is an even function of y. If the <5^°' _a ^(x) are orthogonal to Qq (x) = 
1 then the sum in Eq. (|42p will collapse to its first term and Jo is a constant. It 
is quite a surprise that this actually happens in the interval — 1 < y < 1 for all 
parameter values a for which the integral Jo converges. For example, Io{y) = 
r(a)ir =const. with a rational number r(a) that depends on the exponent a, 
where r(0) = 1/2, r(l) = 3/4, r(2) = 5/8, r(3) = 35/64, r(4) = 3 2 • 7/2 7 , if 
a is a non- negative integer; in general Io(y) = \^7fT(a + 3/2)/T(a + 2). For 
y > |1|, Io(y) varies and deviates from the constant value. From the structure 



(42) 
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of the integral, this anomalous behavior of Jo is rather unexpected. Since 
for a = parity is conserved in the ODE (11 II) . the orthogonality integrals 
O^z/'o are zero an d Or^^p = more generally. Since it is shown in ref. [3] 
that O2J2U f° r ^ ^ v vanish, the Q^~ a \x) polynomials are orthogonal in 
the conventional sense. Since each C/ 0, (x) is proportional to Qf~ a \x), 
the cj; ' a \x) polynomials are orthogonal. This is confirmed by Iq and its 
constancy in the interval — 1 < y < 1 is thus proved. 

The restriction to parameter value a = can be removed: 



-00 a a + 2 [l+y 2 + \^\ a+x 

oo d x e ~acof 1 (x+y y /a) ^ ^ -a. cot" 1 (x+y^/a) 



(43) 



which is neither even in ?/ nor independent of y. Therefore, if we wish to 
find the normalizations of the Q)^~ a ' polynomials we have to split up the 
generating function into its even and odd parts in y and integrate them 
separately, each with the proper power of ct(x) in the orthogonality integral. 



4 Relations with Gegenbauer Polynomials 

The relation of the Romanovski polynomials as complexified Jacobi polyno- 
mials on the unit circle in the complex plane is described in detail in ref. [5]. 
Therefore, we focus here on relations with Gegenbauer polynomials. 

We start with the simplest case of parameter values a = = a, which 
also happens to be relevant for physics [4J, to derive from the generating 
function an expression for Q^° i ~ a '(x) in terms of a finite sum of Gegenbauer 
polynomials. For a = = a, Eq. (j!9p takes the explicit form 

Q(x, y- 0, 0) = (1 + x 2 ) [1 + (x + y(l + x 2 )) 2 ]" 1 = — -^7— 2V (44) 

1 + 2xy + y z {l + arj 

Theorem J^.l. The Q^'°\x) polynomials have the expansion into Gegen- 
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bauer polynomials 

[m/2] 



QT\x) = ml £ (-l)V"CtS(-x), m = 0, 1, . . . . (45) 

n=0 



Proof. For |xy|/|l + 2xy + y 2 \ < 1, the generating function identity 
may be expanded as an absolutely converging geometric series 



-x 2 y 2 ^ 



Q(s, y; 0, 0) = £ v ^ (46) 
Substituting the generating function of Gegenbauer polynomials [7] Cf 1 {x), 



(l-2xy + y 2 )-^ =Y J C\ n+l \x)y\ (47) 

1=0 

we obtain the expansion 

oo oo 

Q(x,y; 0,0) = E(-^V) n E^ (n+1) (-^y 

n=0 i=0 

oo [m/2] 

= £ 2/ m £ (-l)V n C^2(-x), (48) 

m=0 n=0 

where m = / + 2n was used upon interchanging the summations, with [m/2] 
the integer part of m/2. On comparing with Eq. f|T9|) defining the generating 
function Q(y,x; a, —a) we obtain the expansion of the Q^'°\x) polynomials 
as a finite sum of Gegenbauer polynomials of Theorem 4.1. o 

Since Q^°\x) = K m C^(x), this result is also valid for the C£'°\x) 
polynomials. It can be generalized to parameter values a ^ : 

Theorem 4-2. TheQ^' (x) polynomials have the Gegenbauer polynomial 
expansion 

W 2 } / n 1 \ 

n 



Q^- a \ X ) = m £ ( \ ) . 2n crr(-4 (49) 

n=0 



Proof. This relation follows from expanding the generating function 
/ erf a;; \ 1 

Q(x,y;0,-a) = 



1 + x 2 + y 2 a 2 (x) + 2xya(x) J [1 + 2xy + y 2 + x 2 y 2 ] a+1 
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E 



a — 1 \ [x y 



2„,2\n 



n=0 

00 ' -a - 1 



n / [1 + 2xy + y 



21n+a+l 



n,Z=0 x n 

00 [N/2] 



E ( u IW^^V 



—a — 1 



E^E ( n )^cSS hl, (-«) (so) 



AT=0 n=0 

in terms of the binomial series and then again using the generating functions 
of the Gegenbauer polynomials, o 

Theorem 4-3. The Q^' (x) polynomials have the general Gegenbauer 
polynomial expansion 

1 N / N \ \ N I 2 ^ / _ _ 1 \ 

>t ara) W = EL *-iWE a . )^_ + 2 a n +1) (-^)- (5i) 

Proof. Substituting the expansion of Theorem 4.2 into Eq. (1351) in which 
the Gegenbauer polynomials depend only on the parameter a while the 
Q^?' 1 ) (x) depend only on a yields the desired expansion, o 

The Gegenbauer polynomials are well-known generalizations of Legendre 
polynomials. The hyperbolic Gegenbauer ODE 

a(x)y"~(2X + l)xy' + A { l X) y = (52) 

becomes the ODE (11 II) for a = 0, v = I and 2A + 1 = 2(a + u), so the solutions 
of Eq. (1521) are the Qf' 1+1 ' (x) polynomials. In fact, for a = we can 
directly solve the ODE (fTTT) for the Q^~ a \x) polynomial solutions in terms 
of finite power series. 

Theorem 4-4- The Q^'~ a \x) polynomials have the explicit finite power 
series 

Qr .) W = (JT- ^ + 2)^-2^, + 1) 

^ N V ; ^ ^ M 2/i(2a + 2^+1) M 



iV(iV-l) iV(iV-l)---(iV-2 / i + l) 

ai = -7777; — —77, ^ = 



2(2a + 3)' ^ v y 2^!(2a + 3)(2a+l)---(2a + 2/i + l) 

(53) 

N //Vx [iV/2] 

*- a) (^) = E J tt) E (54) 

!/=o v ^ y ^=0 



14 



with from Eq. (1531) and [N/2] denoting the integer part of N/2. 

Proof. Since the proof by mathematical induction is straightforward, we 
just give the results. As the ODE is invariant under the parity transforma- 
tion, x — > —x, we have even and odd solutions. Substituting Eq. fl53l in 
Eq. fl35|) yields the second relation stated in Theorem 4.4. o 

This is also valid for the C^ ' (x) polynomials up to the normalization 

K h 

5 Auxiliary polynomials 

Carrying out the innermost derivative of the Rodrigues formula (pQ), we find 

^(0,-0-1)/ x / , d l ~ l ((?'Wo\ . . 

= aQU ^)-( a +l)-— (— j, (55) 
and are led to define auxiliary polynomials: 

S l+1 (x) = a ^ +1 - ( °'( x )™o( x ) \ / 56 n 
wq(x) dx l \ o~{x) ) 

Example. 

Si(x) = cr'(x), 5*2(2;) = a"a(x) + a'(x)[a — a'(x)(a + 2)], o (57) 

So 

S l (x) = -^-Q\X a - 1 \x)-^. (58) 

Ci — r~ X (jj ~\~ a. 

Applying a derivative to w Si/a l yields 

d l ( a'(x)w Q (x) \ a d /w ( a - a -i) \ I d f w Pi(x) \ 
dx l { a{x) ) ~ a + 1 dx U^- 1 [X) ) a+ldx{ a 1 ) ' l j 

Using the recursive ODEs for Q\l\°~ 1) and P, yields 
tuo(^) cte' \ er(x) / a + 1 

_ gfrr a) (s) 

a + 1 V ; 
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A comparison of Eqs. (I58|I60I) yields 

P l+1 (x) = aa'(x)Q < i a _\ a ' 1 \x) + Qil\ a \x), (61) 

complementing the relation V\{x\ I) = KiC\ a ~ a \x) = Q ( f~ a \x) = Pi(x). 
For Laguerre polynomials <j(x) = x and relation fl55|) corresponds to 

S l+1 (x) = x l+1 e x ^ {^pj = ULT l -\x), (62) 

while Eq. ( |55|) becomes 

lLi(x) = lL[_i(x) — xL i 1 _ 1 (x). (63) 

For Jacobi polynomials <y'(x) = —2x. As — 2x = 1 — x — (1 +x), where 1 ± x 
can be incorporated into the weight functions w(x) = (1 — x) a (l+x) b , there is 
no need to introduce auxiliary polynomials. For example, Eq. fl53|) becomes 

2lP l {a ' b \x) = (a + 0(1 + x)P^f +1 \x) - (b + i)(l - x)PlS 1 ' b \x). (64) 

6 Discussion 

We have used a simple and natural method for constructing polynomials 
Ql?~ a Kx) that are complementary to the C^~ a \x) polynomials and related 
to them by a Rodrigues formula. Similar to the classical orthogonal poly- 
nomials, the Q^~ a \x) appear as solutions of a Sturm-Liouville ordinary 
second-order differential equation and obey Rodrigues formulas themselves. 
On the other hand, and different from the classical polynomials, their infinite 
sets of orthogonality integrals are not the standard ones. These real orthog- 
onal polynomials and their nontrivial orthogonality properties are closely re- 
lated to Romanovski polynomials and to physical phenomena. In summary, 
all basic properties of Romanovski polynomials derive from the Rodrigues 
formula ([1]) except for the orthogonality integrals. 
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